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ABSTRACT In this paper, the design of the control chart when the variable of interest follows the gamma 
distribution under the neutrosophic statistical interval method (NSIM) is proposed. The average run length, 
probability in-control, and probability of out-of-control will be derived using the NSIM. The neutrosophic 
control chart coefficient will be determined by the algorithm under the NSIM. The neutrosophic average run 
length for various shifts and specified parameters will be determined. The efficiency of the proposed control 
chart is discussed using the simulation study and a real example. 


INDEX TERMS Neutrosophic statistics, classical statistics, fuzzy logic, neutrosophic logic, neutrosophic 


average run length. 


I. INTRODUCTION 
High quality is the main aim of well-reputed industries or 
service companies in the modern world. The high-quality 
target can be only achieved by minimizing the non- 
conforming or the defective items through a proper monitor- 
ing of the manufacturing process. The proper monitoring of 
the process can be done using the control charts. The control 
charts analyzed the process, help the industrial engineers, 
and provide information to management in the variety of 
ways such as when the corrective action should be taken, 
unexpected shift in the process and pinpointing the causes of 
variation in the process. Several authors designed a control 
chart for the various reasons, for example, [1] discussed the 
application of the control chart in medicine and epidemiol- 
ogy. Reference [2] applied the control chart in dairy form. 
Detailed applications for healthcare issues can be seen in [3]. 
The Shewhart control charts are designed under the 
assumption that quality of interest follows the normal dis- 
tribution. But, the several manufacturing processes such as 
chemical, semiconductor, cutting tool wear and accelerated 
test the distribution of variable of interest 1s skewed, see 
[4] and [5]. According to [6], the use of the Shewhart control 
chart when data is not collected in the subgroup or when 
the variable of interest follows the skewed distribution may 
mislead the industrial engineers in interpreting the state of the 
control chart. The gamma distribution is an appropriate model 
for the skewed data and used to model the time between 


events. Several authors worked on the control chart charts for 
the gamma distribution, for example, [7] designed the R chart 
for the gamma distribution. Reference [8] designed control 
chart when the quality of interest follows the gamma distri- 
bution. Reference [9] worked on the control chart when the 
time between follows the gamma distribution. Reference [10] 
proposed control for the gamma distribution using the belief 
information. Reference [11] studied the control chart for 
this distribution using a narrow confidence interval method. 
References [12] and [13] designed a chart for the gamma 
distribution under resampling approaches. 

The fuzzy logic is applied in the designing of control charts 
when the proportion of defective parameter is not determined 
value. According to [14] ““The use of a fuzzy approach in 
the design of control charts has allowed to improve the per- 
formance of traditional charts, as well as to be possible of a 
simple approach for the design of control charts for linguistic 
variables with multinomial distribution for both, the univari- 
ate case as for the multivariate case’’. Due to the importance 
of the fuzzy- based control charts in the industry, several 
authors designed control chart based on fuzzy logic includ- 
ing for example, [15] and [16] introduced the fuzzy-based 
algorithm to improve the process. Reference [17] worked on 
fuzzy control chart for skewed distribution. Reference [18] 
discussed the cost model for the Weibull distribution under 
the fuzzy approach. Reference [19] designed control chart 
for the fuzzy random numbers. Reference [20] worked 
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on the fuzzy X-bar chart using non-normal distribution. 
Reference [21] control chart for multinomial control chart 
using fuzzy logic. 

Reference [22] developed the generalized form of fuzzy 
logic is called the neutrosophic logic. The neutrosophic 
logic becomes the traditional fuzzy logic when there are 
indeterminate observations in the sample or in the popu- 
lation. Reference [23] developed the generalized class of 
the classical statistics using the neutrosophic logic is called 
the neutrosophic statistics (NS). The NS reduced to clas- 
sical statistics when all the data is clear and determined. 
References [24] and [25] used the NS to study rock 
measuring issues. References [26]-[28] introduced the NS 
in the area of acceptance sampling plans. Recently, [29] 
introduced the attribute control chart based on the NS. 
Reference [30] proposed the variance control chart using the 
neutrosophic statistical interval method (NISM). More details 
can be seen in [30]. 

The existing control charts for the gamma distribution 
under the classical statistics cannot be applied when the 
observations are fuzzy, indeterminate and unclear. There- 
fore, our aim is to design a control chart to monitor the 
gamma-distributed quality of interest under the uncertainty 
environment. By exploring the literature and the best of our 
knowledge, there is no work on the control chart for the 
gamma distribution under the NISM. In this article, the design 
of the control chart when the variable of interest follows the 
gamma distribution under the NSIM is proposed. The average 
run length, probability in-control and, probability of out-of- 
control will be derived using NSIM. The neutrosophic control 
chart coefficient will be determined the algorithm under the 
NSIM. The neutrosophic average run length (NARL) for 
various shifts and specified parameters will be determined. 
The efficiency of the proposed control chart is discussed 
using the simulation study and a real example. 


ll. NEUTROSOPHIC GAMMA DISTRIBUTION 

Suppose that the neutrosophic life/failure time Ty € 
[TL, Ty], where Tg and Ty present the lower and upper 
failure of indeterminacy interval of an item follows the 
neutrosophic gamma distribution with neutrosophic shape 
parameter ay € [az, ay] and neutrosophic scale parameter 
by € [bL, by]. The neutrosophic probability density func- 
tion (npdf) of the neutrosophic gamma distribution (NGD) is 
given by 





bi’ 1 
f tn) = o aye ONIN: ty, an, by > 0; 


by € [bz, bu] (1) 


where I (x) presents the neutrosophic gamma function, 
see [28]. The corresponding neutrosophic cumulative distri- 
bution function (ncdf) of NGD is given by 


an € [aL, auy], 


iN : 
ay—Le ’N (ty /bNY | 


P(Ty <ty)=1-% > Ty ElTL, Tul, 


j=l i! 
an € [az, au], bn € [bL, byu] (2) 
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Note here that the NGD is the generalization of traditional 
gamma distribution under the classical statistics. The NGD 
becomes the traditional gamma distribution when all the 
observations or the parameters are determined values. The 
neutrosophic mean and neutrosophic variance of NGD are 
given by 


a 
wv = 503 ay € lar, av], by € Ibr, bul (3) 
and 
a 
oy = 53 an € lar, av], by € [bz bu] (4) 
N 


By following [31], the NGD can be transformed to an 
approximately neutrosophic normal distribution when TÈ = 
ii >. Ty € [Tz, Ty]. More details about the neutrosophic 
distributions can be seen in [23] and [26]. The neutrosophic 
mean and neutrosophic variance of Ty, € i. T| are 


given by 


by T(an + 1/3). 


m= Pan) ; ay €[az, ay], by € [bL, bu] 
(5) 
and 
2 
o by T(an+2/3) (bN Tan +1/3)\ 
a ray) ray) | 


an €laz,au], bn €[bt,bu] (6) 
lil. DESIGN OF THE PROPOSED CONTROL CHART 

As mentioned by [31], the transformed vaiable Ty, = 
fi Ty € [TL, Ty] have the symmetry property of the 
neutrosophic normal distribution. The ncdf is valid for integer 
values of ay. We propose the following control chart under 


the NISM when the quality of interest follows the NGD. 


1. Compute Ty = TN ? for the randomly selected item 


from the production process. 

2. Plot Ty, and declare the process out-of-control if Ty, > 
UCLy or Ts < LCLy; where LCLn € [LCL,, LCLy] 
and UCLn € [UCL,, UCLy] are neutrosophic lower 
control limit and neutrosophic upper control limit, 
respectively. 

The proposed control chart under the NISM is the exten- 
sion of the control chart proposed by [8] under the classical 
statistics. The proposed control chart becomes [8] chart when 
no uncertain observations or parameters exist. Suppose now 
that the process is in in-control state at neutrosophic scale 
parameter bon € [boz, boy]. The neutrosophic control limits 
for the in-control process are given by 


3 
bowl (an +4) 
I’ (ay) 


2/3 
bon Fan +2/3) a2 m 
P(ay) ty 


LCLN = UT% = knor» = 


—kyn 


8859 
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bo Paw + 1/3) 
(ay) 
boi T (ay + 2/3) 
+ kn DE Mrs (8) 


[kz, ky] is the neutrosophic control limit 


UCLN = urs + knor: = 


where ky € 
coefficient. 
Let, we define 


LLN 
P(ay+1/3) _ P(ay+2/3) _ (Tes +. Dy 
T(ay) F (ay) (ay) 
ULN 
— | P(ay+1/3) (ay +2/3) _ (hes $ Dy 
7 I(ay ) l (ay) l (ay) 


Therefore, the neutrosophic control limits can be written as 
follows 


LCLy = bj LLN (9) 
UCiy = bl CULN (10) 


Note here that the when the process is shifted only the 
neutrosophic scale parameters of the proposed chart will be 
changed. In practice, the neutrosophic shape parameter is 
usually known and or fixed accordingly. The probability of 
in-control when the process is at bon € [boz, boy ] under the 
NISM is derived as follows 


P?u = P(Ty < LCLylby = bon) 
+P (T$ > UCLy|by = bon) (11) 


or 


—LL? Jaj 
0 O ay—l e N (LLN?) 
Pass- D EEN 


a—le -ULA (UL3, y 


+ o a 


J! 
We will measure the efficiency of the proposed control chart 
under the neutrosophic average run length (NARL) which 


shows on the average when the process is out-of-control is 
defined by 


ARLon = =, ’ 
out, N 
ARLon € [ARLoz,, ARLoy | (13) 


Suppose that due to some special causes of variations, the pro- 
cess has shifted from the targeted bon € [boz, boy] to bin = 
Chon; bin € [biL, biy], where the constant c shows the shift 
in the process. The probability of in-control when the process 
is at bin E [b1L, biy ] under the NISM is derived as follows 
Plan = P (T < LCLn|biw = cbon) 
+P(Ty > UCLy|bin =cbon) (14 
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TABLE 1. The values of NARL when ARLon = 200, 300, 370. 


N [1.9,2.1] 

NARI ns 

4 

3 

2.8 

25 

2.25 [6.23,3.90] 

2 | [7.85,5.08 [9.17,5.79] 
1.9 | [9.37,6.12] [12.079,7.657] 


1.8 [15.048,9.683] 

L7 
1.6 
1.5 


1.4 

1.3 

1.2 [111.95,97.87] | [131.618,115.703] 
11 [224.98,219.282] 

1 [370.011,375.064] 
[461.386,307.94] 

0.75 [404.591,239.851] 

0.7 [340.814,181.489] 
[178.81,79.51] | [224.456,97.608] 

0.5 [109.06,39.69] 

0.4 

0.3 

[22.609,5.004] 
[4.35,1.41] [5.94,1.58] 
[2.22,1.06] 





or 


-UL}, UL?, j 


-1e c © (=i 
ae — s 


J! 
(15) 








aaa ) 


ay —1 c 
Plen = 1- pe 


The NARL for the shifted process is defined as 


ARL in = ARL yn € [ARL1,,ARL iy] (16) 


N 
Suppose that roy be the specified value of ARLon. The 
values of ARL yn for various shift c, ay and by are dis- 
played in Tables 1-2. From Tables 1-2, we note that for the 
fixed values of ay € [3,5], by € [1.9,2.1] and c, the 
indeterminacy interval range in ARL1n increases as ARLon 
increases. We also note that indeterminacy interval range in 
ARL yn increases when ay € [3,5]; by € [1.9, 2.1] change 
to ay E€ [5, 10]; by € [1.45, 1.55]. 

The following neutrosophic algorithm is used to find 
kn € [kz, ky] and ARLyn € [ARL1,, ARLiy]. 


1. Specify the values of ay € [az,ay],bn € [bL, bu] 
and c. 
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TABLE 2. The values of NARL when ARLon = 200, 300, 370]. 


[2.77,2.79] [2.89,2.92] [2.96,2.98] 


ae tid) —_f1e) a 


[1.45,1.55] [1.45,1.55] [1.45,1.55] 


NARL 
[1.33,1.04] [1.40,1.05] 


[1.37,1.05] 


B a a 
2.8] (2.101.295) | 2.27134] | 2.36137] 
2.5 | [2.66152] | [2.92,1.60] | [3.06,1.64] | 
2.25 | [3.481.871] | [3.89,2.01] | 4.13,2.08] | 

Aa a e 

18 [748,378] o oe 
i a locas Teen 
16 | 12.74,6.59] | [15.57,7.79] | [17.28,8.49] | 
1.5 | [17819.55]  |[22.23,11.57] | [24.94,12.75] | 
14 | [26.39,15.03] | [83.78,18.74] | [38.38,20.94] | 
13 | [41.83,26.22] | [55.18,33.81] | [63.67,38.42] | 
1.2 | [71135157]  |[97.38,69.31] | [114.49,80.41] | 
14 | [125.6,111.45] | [179.98,157.93] | [216.67,188.06] | 
1 | [200.04,203.21] | [300.02,302.98] | [370.01,370.18] _| 
0.8 | [163.05,92.11] | [245.46,132.99] | [303.72,160.07] | 
0.75 | [128.44,61.92] | [191.91,88.21] | [236.61,105.49] | 
0.7 | [98.27,40.97] | [145.75,57.52] | [179.08,68.31] | 

0.5 | [27.61,7.18] [39.46,9.302] [47.63,10.62] 
0.4 | [12.88,3.04] | [17.84,3.69] | [21.22,4.09] | 
03) (5.43148)  |[7.17,1.64] | [8321.75] | 
0.25 | (3.41116) | [4341.23] | [4.961.227] | 
0.15 | [1.40,1.00] | [1.58,1.00] | [1701.00] | 
04 [[1.059,1.00]  |[1.10,1.00]  |[1.13,1.00] | 
0.05 | [1.00,1.00] | [1.00,1.00] | [1.00,1.00] | 














2. Determine kyn € [kL, ky] such that ARLon => ron. 

3. Several combinations of ky € [kz, ky | will exist where 
ARLon = ron. 

4. Choose that values of kn € [kL, ky] where ARLon is 
ron Or very close to it. 

5. Use the selected kn € [kz, ky] to find ARLin € 
[ARL1,, ARL1y | at various values of c. 


IV. ADVANTAGES OF THE PROPOSED CONTROL CHART 
The efficiency of the proposed control chart under NISM 
with the control chart under the classical statistics will be 
discussed in this section. For the fair comparison, we will 
fix the same values of specified parameters for both control 
charts. 


A. THEORETICAL COMPARISON 

In this section, we show the efficiency of the proposed control 
chart under NISM over the existing control chart proposed 
by [8] under the classical statistics. Both control chart will 
be compared in NARL. The control chart having the smaller 
values of NARL is known as the most efficient chart. The 
smaller the NARL, quicker the indication about the shift 
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TABLE 3. Comparison of proposed chart and existing chart. 




















ARLs for ARLin of the | ARLs for 
ARLin of the existing proposed existing 
c proposed chart chart c chart chart 
sF paz 
3 [3.2,2.07] [224.98,219.28] 230.25 
2.8 [3.71,2.37] [370.01,375.06] , 
2:5 [4.92,3.08] [461.39,307.94] 
225 [6.67,4.15] [404.59,239.85] 
2 [9.94,6.24] [340.81,181.49] 
1.9 [12.08,7.66] [224.46,97.61] 
1.8 [15.05,9.68] [136.35,48.22] 
i7 [19.32,12.7] [74.96,21.46] 





1.6 [25.69,17.4] 
1.5 | [85.64,25.13] [22.61,5] 


1.4 [51.98,38.7] 52.83 | 0. L [7.05,1.71] 7.20 
1.3 | [80.27,64.26] 81.74 [3.25,1.13] 3.30 


in the process. Also, as pointed out by [24] that a method 
which has provided the interval range rather than the deter- 
mined value of the parameters is known as most effective and 
adequate under the uncertainty environment. The values of 
ARLyn € [ARL1,, ARL jy | of the proposed control chart and 
the existing control chart when ay € [3,5], by € [1.9, 2.1] 
and ARLon € [370, 370] are shown in Table 3. 

From Table 3, it can be observed that the proposed control 
chart has smaller values of NARL as compared to the chart 
proposed by [8] at all levels of c. For an example, when 
c = 0.7, the indeterminacy interval of NARL is ARLin € 
[340.81, 181.49] while the existing control chart has a deter- 
mined value of ARL is 350.95. Under the indeterminate 
environment, the proposed control chart indicates that the 
process will be shifted between 181% sample and 340" sam- 
ple. On the other hand, the existing chart indicates that the 
process will be shifted at the 350" sample. By comparing 
both control chart, it is concluded that the proposed control 
chart not only provides the quick indication about the shift in 
the process but also provide the indeterminacy interval which 
is required under the uncertainty environment. Therefore, 
the theory of the proposed control chart matches with [24]. 


[35.58,8.41] 




















B. THEORETICAL COMPARISON 

Now, we will compare the proposed control chart under the 
NISM with the existing control chart using the simulated 
data. The data is generated from the NGD when ay € 
[3,5], by € [1.9, 2.1] and ARLon € [370, 370]. First 20, 
observations are generated from the NGD when the process 
is in-control state and next 30 observations are generated from 
the NGD at c = 1.5. At these specified parameters, the tabu- 
lated ARL yn E€ [35.642, 25.132] which means, the first out- 
of-control can be expected between 25" and 35" samples. 
We determined statistic T = ite >. Ty € [TL, Ty] and plot 
it on control chart in Figure 1. From Figure 1, it can be seen 
that the proposed control chart under NISM detects shift at the 
25'" sample. We also plotted the Ty, under the classical statis- 
tics on the control chart in Figure 2. From Figure 2, we note 
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= 
UCLU=3.1076 
La 
™ 
« =. 
LCLL=0./174 
co 





0 10 20 30 40 50 
sample Number 


FIGURE 1. The proposed control chart for simulated data. 


~I 
Co 
UCLL=2 5437 
Cl 
“KE 
LCL=0.4057 
I 





0 10 20 30 40 50 


Sample Number 


FIGURE 2. The existing control chart for the simulated data. 


that the existing control chart does not provide any signal 
about the shift in the process. By comparing Figure 1 with 
Figure 2, we conclude that the proposed control chart is more 
efficient in detecting a shift in the process. 


V. CASE STUDY 

In this section, we will discuss the application of the pro- 
posed control chart under NISM in the healthcare depart- 
ment. The large hospital in Saudi Arabia wants to apply 
the proposed control chart for the monitoring of the urinary 
tract infections (UTIs). The hospital wants to monitor the 
increase or decrease infection rate in those patients who 
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TABLE 4. The neutrosophic ATis data. 
































Sr# Ty Tw Sr# Ty 
1 [6.53,12.47] 
2 | [3.57,15.55] _| [1.52,2.49] | 27 | [6.85,12.13] 
3 [8.08,22.69] 
4 
5 
6 | [11.45,13.18] | [2.25,2.36] | 31 | [6.81,17.25] | [1.89,2.58] 
7 | [9.20,15.18] | [2.09,2.47] | 32 | [4.42,12.53] | [1.64,2.32] 
8 | [5.51,9.77] | [1.76,2.13] | 33 | [6.53,13.96] | [1.86,2.40] 
9 
10 
11 | [7.35,14.89] | [1.94,2.46] 
12 | [5.62,11.09] | [1.77,2.23] 
13 | [8.38,16.72] | [2.03,2.55] 
15 | [4.90,23.67] | [1.69,2.87] | 40 | [5.74,11.95] | [1.79,2.28] 
16 | [4.45,14.68] | [1.64,2.44] 
17 | [7.11,16.44] | [1.92,2.54] [5.10,10.10] 
18 | [9.37,15.95] | [2.10,2.51] [9.11,24.54] 
19 | [12.00,16.38] | [2.28,2.53] 
20 | [7.41,16.62] | [1.95,2.55] 
21 | [10.64,15.15] | [2.19,2.47] | 46 | [7.9,11.43 [1.99,2.25] 
22 | [6.63,11.21] | [1.87,2.23] | 47 | [3.62,13.00] | [1.53,2.35] 
23 | [2.87,14.27] | [1.42,2.42] | 48 | [5.01,13.62] | [1.71,2.38] 
24 | [6.87,10.37] | [1.90,2.18] 
25 | [6.16,18.85] | [1.83,2.66] 





had acquired UTIs during their stay in the hospital. The 
microbes which are too small are the main reason for UTIs 
infection. As pointed out by [32] that the uncertainty exists 
while collecting the UTIs data based on the clinical diagno- 
sis. Therefore, under the uncertainty environment, it is not 
possible that all the UTIs data is determined. So, for the UTIs 
data having some uncertain observations, the existing control 
chart using under the classical statistics cannot be applied for 
the monitoring of infection rate adequately and effectively. 
By following, [6], the data is collected in the number of days 


which follows the NGD with ay € [5, 10], by € [1.45, 1.55]. 


The data along with the statistic T = ii >. Ty € (TL, Tuy] 


is Shown in Table 4. 

Let ARLon € [370, 370], an € [5, 10], bon € [1.45, 1.55] 
and ky € [2.96, 2.981]. The neutrosophic control limits using 
this information is computed as follows 


1/3 
boa (an ++ 5) 


Loree 
I (an) 
2/3 2 
k bonT (av + 3) 2, 
NV FE (ay) OTR 
LCLy € [1.0489, 2.7361] 
1/3 
ban, (ay +1/3 
UCL = ON (ay / ) 
I (ay) 
2/3 
sp ky, [PONTON a 
[(ay) ty 


UCLN € [1.6863, 2.7361] 
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= 
=I 
LO 
ii UCLU=3.2449 
= 
w] 
LD 
(ih 
« i= 
O 
C 
LO 
= 
= LCLL=1.0489 
Wo 
0 10 20 30 40 50 
Sample Number 


FIGURE 3. The proposed control chart for UTIs data. 


35 40 


3.0 


UCLL=2. 7456 


15 20 25 


1.0 


0.5 





0 10 20 30 40 30 


Sample Number 


FIGURE 4. The existing control chart for UTIs data. 


The values of Tx, are plotted on the control limit in Figure 3. 
From Figure 3, it can be read that several values of Ty lie in 
the indeterminacy interval of neutrosophic lower and upper 
control limits. We also note that sample 11" and 21° are close 
to LCL, which needs the hospital management attention. 
On the other hand, the control chart in Figure 4 under the 
classical statistics does not show any point indeterminacy 
interval. 


VI. CONCLUDING REMARKS 

We presented a control chart when the quality of interest 
follows the NGD. The proposed control chart is the gener- 
alization of the control chart under the classical statistics. 
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The proposed control chart can be applied when the 
data or the parameters are indeterminate values rather than the 
exact or determined values. The efficiency of the proposed 
control chart is discussed using the simulation and the real 
example. By compression of the proposed chart with the 
existing control chart, we conclude that the proposed con- 
trol chart performs better under the uncertainty environment. 
We recommend that the industrial engineers and health per- 
sonals should apply the proposed control chart for monitoring 
the data obtained from the sample or the population having 
incomplete information. The proposed control chart for the 
marine big data and ocean big data can be considered as future 
research. 


ACKNOWLEDGEMENTS 

The authors would like to thank the editor and the reviewers 
for their valuable suggestions to improve the quality of this 
manuscript. They would like to thank the Deanship of Scien- 
tific Research for the technical support. 


REFERENCES 


[1] L. H. Sego, “Applications of control charts in medicine and epidemiol- 
ogy,” Ph.D. dissertation, Virginia Tech, Blacksburg, VA, USA, 2006. 

[2] B. Stewart, R. James, K. Knowlton, M. McGilliard, and M. Hanigan, 
“An example of application of process control charts to feed management 
on dairy farms,” Prof. Animal Sci., vol. 27, no. 6, pp. 571-573, 2011. 

[3] G. Suman and D. Prajapati, “Control chart applications in healthcare: 
A literature review,” Int. J. Metrol. Qual. Eng., vol. 9, May 2018, Art. no. 5. 

[4] D.-S. Bai and I. Choi, “X and R control charts for skewed populations,” 
J. Qual. Technol., vol. 27, no. 2, pp. 120-131, 1995. 

[5] K. Derya and H. Canan, “Control charts for skewed distributions: Weibull, 
gamma, and lognormal,” Metodol. Zvezki, vol. 9, no. 2, p. 95, 2012. 

[6] E. Santiago and J. Smith, “Control charts based on the exponential dis- 
tribution: Adapting runs rules for the ¢ chart,’ Qual. Eng., vol. 25, no. 2, 
pp. 85-96, 2013. 

[7] I.M. Gonzalez and E. Viles, “Design of R control chart assuming a gamma 
distribution,” Econ. Qual. Control, vol. 16, no. 2, pp. 199—204, 2001. 

[8] S. H. Sheu and T. C. Lin, “The generally weighted moving average control 
chart for detecting small shifts in the process mean,” Qual. Eng., vol. 16, 
no. 2, pp. 209-231, 2003. 

[9] C. W. Zhang, M. Xie, J. Y. Liu, and T. N. Goh, “A control chart for the 
Gamma distribution as a model of time between events,” Int. J. Prod. Res., 
vol. 45, pp. 5649-5666, Mar. 2007. 

[10] M. Aslam, N. Khan, and C.-H. Jun, “A control chart using belief informa- 
tion for a gamma distribution,’ Oper. Res. Decis., vol. 26, no. 4, pp. 5-19, 
2016. 

[11] S. Hao, S. Huang, and J. Yang, “Design of Gamma control charts based 
on the narrowest confidence interval,” in Proc. IEEE Int. Conf. Ind. Eng. 
Eng. Manage. (IEEM), Dec. 2016, pp. 219-223. 

[12] N. Khan, M. Aslam, L. Ahma, and C.-H. Jun, “A control chart for gamma 
distributed variables using repetitive sampling scheme,” Pakistan J. Stat. 
Oper. Res., vol. 13, no. 1, pp. 47-61, 2017. 

[13] M. Aslam, O.-H. Arif, and C.-H. Jun, “A control chart for gamma distri- 
bution using multiple dependent state sampling,” Ind. Eng. Manage. Syst., 
vol. 16, no. 1, pp. 109-117, 2017. 

[14] M. N. P. Fernandez, “Fuzzy theory and quality control charts,” in Proc. 
IEEE Int. Conf. Fuzzy Syst. (FUZZ-IEEE), Jul. 2017, pp. 1-6. 

[15] H. Rowlands and L. R. Wang, “An approach of fuzzy logic evaluation and 
control in SPC,” Qual. Rel. Eng. Int., vol. 16, no. 2, pp. 91—98, 2000. 

[16] S. M. El-Shal and A. S. Morris, “A fuzzy rule-based algorithm to improve 
the performance of statistical process control in quality systems,” J. Intell. 
Fuzzy Syst., vol. 9, pp. 207—223, Jan. 2000. 

[17] R.Intaramo and A. Pongpullponsak, “Development of fuzzy extreme value 
theory control charts using a-cuts for skewed populations,’ Appl. Math. 
Sci., vol. 6, no. 117, pp. 5811-5834, 2012. 

[18] P. Charongrattanasakul and A. Pongpullponsak, “Economic model for 
fuzzy Weibull distribution,” in Proc. Int. Conf. Appl. Statist., 2014, p. 155. 


8863 


IEEE Access 


[19] S. Mojtaba Zabihinpour, M. Ariffin, S. H. Tang, and A. Azfanizam, 
“Construction of fuzzy X-S control charts with an unbiased estimation of 
standard deviation for a triangular fuzzy random variable,” J. Intell. Fuzzy 
Syst., vol. 28, no. 6, pp. 2735-2747, 2015. 

[20] C. Panthong and A. Pongpullponsak, “‘Non-normality and the fuzzy theory 
for variable parameters control charts,” Thai J. Math., vol. 14, no. 1, 
pp. 203-213, 2016. 

[21] M. Mashuri and M. Ahsan, “Perfomance fuzzy multinomial control chart,” 
J. Phys., Conf. Ser., vol. 1028, no. 1, p. 012120, 2018. 

[22] F. Smarandache, “‘Neutrosophic logic—A generalization of the intuition- 
istic fuzzy logic,” in Multispace & Multistructure. Neutrosophic Transdis- 
ciplinarity (100 Collected Papers Science), vol. 4. Hanko, Finland: North- 
European Scientific, 2010, p. 396. 

[23] F. Smarandache. (2014). Introduction to Neutrosophic Statistics: Infinite 
Study. [Online]. Available: https://arxiv.org/pdf/1406.2000 

[24] J. Chen, J. Ye, and S. Du, “Scale effect and anisotropy analyzed for neutro- 
sophic numbers of rock joint roughness coefficient based on neutrosophic 
statistics,” Symmetry, vol. 9, no. 10, p. 208, 2017. 

[25] J. Chen, J. Ye, S. Du, and R. Yong, “Expressions of rock joint roughness 
coefficient using neutrosophic interval statistical numbers,’ Symmetry, 
vol. 9, no. 7, p. 123, 2017. 

[26] M. Aslam, “A new sampling plan using neutrosophic process loss consid- 
eration,’ Symmetry, vol. 10, no. 5, p. 132, 2018. 

[27] M. Aslam and M. A. Raza, “Design of new sampling plans for multiple 
manufacturing lines under uncertainty,’ Int. J. Fuzzy Syst., pp. 1-15, 
Oct. 2018, doi: 10.1007/s408 15-018-0560-x. 

[28] M. Aslam and O. Arif, “Testing of grouped product for the weibull 
distribution using neutrosophic statistics,’ Symmetry, vol. 10, no. 9, p. 403, 
2018. 

[29] M. Aslam, R. A. R. Bantan, and N. Khan, “Design of a new attribute 
control chart under neutrosophic statistics,” Int. J. Fuzzy Syst., pp. 1-8, 
Nov. 2018, doi: 10.1007/s40815-018-0577-1. 

[30] M. Aslam, N. Khan, and M. Khan, “Monitoring the variability in the 
process using neutrosophic statistical interval method,” Symmetry, vol. 10, 
no. 11, p. 562, 2018. 

[31] E. B. Wilson and M. M. Hilferty, ““The distribution of chi-square,” Proc. 
Nat. Acad. Sci. USA, vol. 17, no. 12, pp. 684—688, 1931. 

[32] R.A. Taylor, C. L. Moore, K.-H. Cheung, and C. Brandt, “Predicting uri- 
nary tract infections in the emergency department with machine learning,” 
PLoS ONE, vol. 13, no. 3, p. e0194085, 2018. 


es 


Wi o m 






MUHAMMAD ASLAM received the M.Sc. and 
M.Phil. degrees in statistics from Government 
College University (GC University) (Lahore), 
in 2004 and 2006, respectively, and the Ph.D. 
degree in statistics from the National Col- 
lege of Business Administration and Economics, 
Lahore, in 2010, under the kind supervision of 
Dr. M. Ahmad. He was a Lecturer of statistics 
with the Edge College System International, from 
2003 to 2006. He was a Research Assistant with 
the Department of Statistics, GC University (Lahore), from 2006 to 2008. 

In 2009, he joined Forman Christian College as a Lecturer, where he was an 
Assistant Professor, from 2010 to 2012. He was an Associate Professor with 
Forman Christian College, from 2012 to 2014. He was an Associate Professor 
of statistics with the Department of Statistics, Faculty of Science, King 
Abdulaziz University, Jeddah, Saudi Arabia, from 2014 to 2017. He taught 
summer course as a Visiting Faculty Member of statistics at Beijing Jiaotong 
University, China, in 2016. He is currently a Full Professor of statistics with 
the Department of Statistics, King Abdulaziz University. He is appointed as 
an External Examiner for 2016/2017—2018/2019 triennium at The University 
of Dodoma, Tanzania. He has published more than 300 research papers in 
national and international well-reputed journals, including the IEEE Access, 






8864 


M. Aslam et al.: Design of a Control Chart for Gamma Distributed Variables Under the Indeterminate Environment 


the Journal of Applied Statistics, European Journal of Operation Research, 
Information Sciences, the Journal of Process Control, the Journal of the 
Operational Research Society, Applied Mathematical Modeling, the Inter- 
national Journal of Advanced Manufacturer Technology, Communications 
in Statistics, the Journal of Testing and Evaluation, and the Pakistan Journal 
of Statistics. His papers have been cited more than 2200 times with an 
h-index of 25 and an 1-10 index of 68 (Google Coalitions). His papers 
have been cited more than 1000 times with an h-index of 19 (Web of 
Science Coalitions). He has authored one book published in Germany. He has 
been an HEC approved Ph.D. Supervisor, since 2011. He has supervised five 
Ph.D. theses, more than 25 M.Phil. theses, and three M.Sc. theses. He is 
supervising one Ph.D. thesis and more than five M.Phil. theses in statistics. 
His current research interests include reliability, decision trees, industrial 
Statistics, acceptance sampling, rank set sampling, neutrosophic statistics, 
and applied statistics. He is amember of the Editorial Board of the Electronic 
Journal of Applied Statistical Analysis, the Asian Journal of Applied Science 
and Technology, and the Pakistan Journal of Commence and Social sciences. 
He is also member of the Islamic Countries Society of Statistical Sciences. 
He received the Chief Minister’s Punjab Merit Scholarship for his M.Sc. 
degree, the Governor’s Punjab Merit Scholarship for his M.Phil. degree, 
the Meritorious Services Award in research by the National College of Busi- 
ness Administration and Economics, Lahore, in 2011, the Research Produc- 
tivity Award by the Pakistan Council for Science and Technology, in 2012, 
the King Abdulaziz University Excellence Award in scientific research for 
the paper entitled "Aslam, M., Azam, M., Khan, N. and Jun, C.-H. (2015). 
A New Mixed Control Chart to Monitor the Process, International Journal 
of Production Research, 53 (15), 4684—4693, and the King Abdulaziz Uni- 
versity Citation Award for the paper entitled "Azam, M., Aslam, M. and 
Jun, C.-H. (2015). Designing of a hybrid exponentially weighted moving 
average control chart using repetitive sampling, International Journal of 
Advanced Manufacturing Technology, 77:1927-1933, in 2018. His name 
listed in second position among statistician at the Directory of Productivity 
Scientists of Pakistan, in 2013. His name listed in first position among 
statistician at the Directory of Productivity Scientists of Pakistan, in 2014. 
He got 371 positions in the list of top 2210 profiles of the Scientist of Saudi 
Institutions, in 2016. He is selected for the Innovative Academic Research 
and Dedicated Faculty Award 2017 by SPE, Malaysia. He is a Reviewer of 
more than 50 well-reputed international journals. He has reviewed more than 
140 research papers for various well-reputed international journals. 


RASHAD A. R. BANTAN received the Ph.D. degree from the Geology 
Department, Royal Holoway, University of London, in 1999. He is currently 
an Associate Professor with the Department of Marine Geology, Faculty of 
Marine Sciences, King Abdulaziz University. He published several papers 
in the well-reputed journals. He is interested to apply statistical methods in 
marine technology. 


NASRULLAH KHAN received the Ph.D. degree in statistics from the 
National College of Business Administration and Economics, Lahore, 
in 2014. He is currently an Assistant Professor in statistics with the Depart- 
ment of Social Science, University of Veterinary and Animal Sciences, 
Jhang. He has published his research work in the field of control chart and 
sampling plans. His research interest includes industrial statistics. 


VOLUME 7, 2019 


